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^ , ABSTRACT 

In this first paper of a series on the structure of boxy and peanut-shaped (B/PS) bulges, 
' itn-band observations of a sample of 30 edge-on spiral galaxies are described and dis- 

cussed, i&i-band observations best trace the dominant luminous galactic mass and 
are minimally affected by dust. Images, unsharp-masked images, as well as major-axis 
^ I and vertically-summed surface brightness profiles are presented and discussed. Galax- 

. ies with a B/PS bulge tend to have a more complex morphology than galaxies with 

other bulge types, more often showing centered or off-centered X structures, secondary 
maxima along the major-axis and spiral-like structures. While probably not uniquely 
related to bars, those features are observed in three-dimensional N-body simulations of 
barred discs and may trace the main bar orbit families. The surface brightness profiles 
, of galaxies with a B/PS bulge are also more complex, typically containing 3 or more 

clearly separated regions, including a shallow or flat intermediate region (Freeman 
Type II profiles). The breaks in the profiles offer evidence for bar-driven transfer of 
angular momentum and radial redistribution of material. The profiles further suggest a 
Q ' rapid variation of the scaleheight of the disc material, contrary to conventional wisdom 

^ I but again as expected from the vertical resonances and instabilities present in barred 

C/3 i discs. Interestingly, the steep inner region of the surface brightness profiles is often 

' shorter than the isophotally thick part of the galaxies, itself always shorter than the 

flat intermediate region of the profiles. The steep inner region is also much more promi- 
nent along the major-axis than in the vertically-summed profiles. Similarly to other 
recent work but contrary to the standard 'bulge -I- disc' model (where the bulge is both 
thick and steep), we thus propose that galaxies with a B/PS bulge are composed of a 
thin concentrated disc (a disc-like bulge) contained within a partially thick bar (the 
B/PS bulge), itself contained within a thin outer disc. The inner disc likely formed 
secularly through bar-driven processes and is responsible for the steep inner region of 
the surface brightness profiles, traditionally associated with a classic bulge, while the 
bar is responsible for the flat intermediate region of the surface brightness profiles and 
the thick complex morphological structures observed. Those components are strongly 
coupled dynamically and are formed mostly of the same (disc) material, shaped by the 
weak but relentless action of the bar resonances. Any competing formation scenario 
for galaxies with a B/PS bulge, which represent at least 45 per cent of the local disc 
galaxy population, must explain equally well and self-consistently the above morpho- 
logical and photometric properties, the complex gas and stellar kinematics observed, 
and the correlations between them. 
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1 INTRODUCTION 

Since the work of iDavies et alJ 1^8^, bulges liave largely 
been considered as low-luminosity elliptical galaxies, sug- 
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gesting a rapid formation dominated either by merg- 
ers a n d a ccretion of e xternal mat e rial ( e.g. iKauffmannI 
| l996l: iKaiif fmann. Cha riot fc White! 11999: but se e also 
bekki fc Chibal I2OOII: lAguerri. Balcells fc Peletieil H 
or possibly by dissipative g r avitat i onal collapse (e 
Eggen. Lynden-Bell fc Sandaed Il962l: ICarlberd Il984alf 



Steinmetz fc Miilledll99,4 but see also lNoguchilll99?j . Il99ty 



Over the last decade, however, there has been much crit- 
icism of this idea. In particular, the realization that most 
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Alternative models where bulges grow secularly (i.e. 
over a long timescale and in relative isolation) have 
also been developed and st udied, many of thom bar- 
drive n (e.g.jPfen niger fc Norman^l990j:jl ^iedli fc Benz-1993. . 
1199,4 iNorman. Sellwood fc HasanI Il99fih. and much ob- 
servational data sugport them fe.g. IWozniak et alJ Il995l : 
lErwin fc Sp arkc 200^. 

Of primary concern here, several pieces of evidence 
point to the identification of most boxy and peanut-shaped 
(B/PS) bulges in edge-on spiral galaxies with the bars of 
barred spirals. In A'^-body simulations, whenever a disc 
galaxy forms a bar, a B/PS bar/b ulge develops soon 
after. This process was studied firs t bv [ Combes ^^^mde^ " 



) , later on by ICombes et al ] Jl990ir and Raha et al 
1991 ), and more rece n tly b y jMihos et alj ^ lll99 "" 



^thanassoula fc Misirioti;^ ("2002'), 'Athanassoula ('200: 
1 2003. 2005) and Martinez- Valoucsta, Shlosman fc HcUerj 

The observed incidence of B/PS bulges is consistent 
with that expected if they a re associated with relatively 
strong bars. Recent work by iLutticke. Dettmar fc PohlenI 
i2000allbl) demonstrates that 45 per cent of all bulges 
are B/PS, while amongst those the exact shape of the 
bulge depends mainly on the viewing angle to the bar. 
As shown by the numerical simulations, true peanuts 
are bars seen side-on, i.e. with the major-axis of the 
bar roughly perpendicular to the line-of-sight. For less 
favourable viewing angles, the bulge/bar looks boxy, 
and if the bar is seen end-on it looks almost spheri- 
cal. Stronger bars also lead to more pro minent peanut 
shapes , as demonstrated observat ionally (e.g. |lAitticke_et_aI 
20001 ) and theoretically (e.g. [ Athanassoula fc Misirioti^ 



200i-]Bureau fc Athanassoulall20oi) . 

The kinematics of discs harboring a B/PS bulge, 
as measured from both ionized-gas emission lines and 
stellar absorption lines, show the behaviour expected of 
barred spirals viewed edge- on. This has been demon - 
strated by the observation s of iKu iiken fc Merrifiel d| lll995l) , 
iMerrifiel fc Kui ikcn (1999). 'B ureau fc FreemanI lll997lTT999l) 
and ICh un^T^ Bureau (2004|[) . an d by the modeling of 
lAthanassoula fc Bureaul l|l999t) and lBureau fc Athanassoiiial 
lll999j .l200^~ Similar tests are now also available for face-on 
bars (Dcbattista ct al. 2005). 

There is thus ample evidence that edge-on galaxies with 
a B/PS bulge are simply barred disc galaxies, and that 
the B/PS bulge s themselves repres ent the thickest parts 
of the bars (see lAthanassoulal 120051 for a review of all ar- 
guments). Yet there has been little convincing evidence 



for this fro m surface photometry alone, the best work be- 
ing that of iLutticke et al ] J2000allbl) . Early studies of non- 
spheroidal bulges used mostly optical images (e.g. Ijarvisi 
ll987l:IShaw fc Gilmorelll98gl : IShawli993) and the interpreta- 
tion was often hampered by the large amount of extinction. 
As we will show in this paper, it truly takes the combination 
of N-body simulations and orbit studies with near-infrared 
(NIR) images to derive direct and convincing photometric 
evi dence relating B/PS bulges and bars. 

[Skokos^atsi^^^Mhanag soula l (|2OO2a|0) and 
iPatsis. Skokos fc Athanassoulal 120021 bOOSalT studied 
the orbital structure of three-dimensional (3D) ba rs ex- 
haust ively (but see also j^fcnnigcr 1984; Pfcnnigor fc Friedli 
I1991I) . They find families of orbits which can not only 
provide the backbone of the boxy and peanut shapes, but 
can also cause local enhancements within the disc itself. 
Since extinction is far less important in the NIR, /('-band 
images are the ideal tool to study the morphology of 
galaxies with a B/PS bulge, to look for similarities with 
barred orbital structures. 

In this paper, the first of a series, we study a sam- 
ple of 30 edge-on spiral galaxies, most of them with a 
B/PS bulge, for which we have obtained high-quality K- 
band im ages. Much complementary data exist for this sam- 
ple (e.g. iBureau fc FreemanM l997i Il999l: IChung fc Bureaul 
l2004l : lBureau fc C hung 2006). but the primary goal here is to 
study the morph ology of the B/PS struc tures and their host 
discs, similarly to lLiitticke et all i2000bl) . As we shall see, we 
find several features ressembling closely those expected of 3D 
bars: bulges with X shapes, secondary disc enhancements, 
inner rings, etc. A second paper discusses scaleheight varia- 
tions in the sam e galaxies iAthanassoula. A ronica fc Bureaul 
j200j, hereafter iPaoer lit but see also lATonic^^^im200^ 

We present our sample in § (21 discuss the observations 
and data reduction in §0 and then describe and discuss the 
resulting images and surface brightness profiles in § 21 and 
1^1 respectively. We examine the direct consequences of our 
results in §|^and conclude briefly in §[7| 



2 SAMPLE DESCRIPTION 



We study here the s a mple of iBureau fc FreemanI il999 
and IChung fc Bureaul (jgOOi), which was drawn from the 
catalo g s of g alaxie s with a B/PS bulge o f Ijaryij lll98d) . 
IShawl l|l987l) and Ide Souza fc dos Aniod il987r) . as well 
as from galaxies with extreme axial ratios (a /b 7; 
iKarachentsev. Karachentseva fc Parnovskvl Il993l) . We in- 
clude galaxies with extended as well as no or confined 
ionized-gas emission, even though the latter were not specifi- 
cally studied by Bureau fc Freeman ( 1999). The sample thus 
consists of 30 edge-on spirals, all of which are visible from 
the south (5 < 15°) and have a bulge larger than Q'.Q in 
diameter (to be able to identify bars kinematically in mod- 
erate seeing). 24 of t he bulges wore classified as having a 
B/PS bulge by IBureau fc Freeman (1999) and 6 constitute 
a 'control' sample and have varied morphologies. 

Basic galaxy properties are listed in Table [H along 
with a classification of the bulges' shape from lLiitticke et al.l 
i2000Ej) . Although still relying on a visual inspection 
of optical contour plots, Liitticko ct al. (2000a) homoge- 
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Table 1. Galaxy sample. 



Galaxy R.A.(2000) Dec.(2000) Type Bulge M'i^ D25 Vq 

J, ° ' " _i 

mag mag arcmin km s 
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Notes: All parameters are from LEDA (Lyon-Meudon Extrag alactic Database) , except the morph ological type JjarvisI Il986l : 
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B/PS bulges and bulge type 4 spheroidal bulges. ESO 151-G004's redshift and corrected absolute B magnitude are from lChung fc Bureaul 
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neously classified over 1350 edge-on galaxies, including 
all our objects. Only 4 bulg es w ere classified d i fferentl y 
by iBureau fc FreemanI Jl999l) and iLiitticke et"ai] ||200o3). 
A complete comparison of the classifications as well as 
of the ionized-gas and stellar kinematics is p rovided in 
IChung fc Bureaul J2004l) . iBureau fc Chun3 ll2006h discuss in 
more depth 3 objec ts where counter-ro tating ionized-gas was 
discovered (see also Chung ct al.""2006'). As demonstrated by 
ILiitticke et al... (.2000a) . the classification of the bulges' shape 
is robust between the optical and NIR, but NIR images are 
essential to study their detailed structure, due to the fre- 
quent and significant extinction from dust in the equatorial 
plane ('e.g. ILiitticke et al? ' 2000b ). All sample galaxies there- 
fore also have D25 ^ 7' (where D25 is the diameter at the 
25 mag arcsec"^ isophotal level in B), which allows for the 
NIR imaging discussed here to be acquired in a reasonable 
amount of time with a small-field camera. 



3 OBSERVATIONS AND DATA REDUCTION 

The observations were carried out at the 2.3m telescope 
of Siding Spring Observatory over 8 runs from 1996 Jan- 
uary to 1997 March, for a total of 24 nights. Mounted at 
the //18 Cassegrain foc us, the CASPIR ^ (Cryogenic Ar- 
ray Spectrometer /Imager: [McGregor et alll994l) instrument 
was used in direct imaging mode with a 0.5 arcsec pix~^ 
scale. Its 256 x 256 pix^ SBRC InSb CCD (30^m pixels) 
then yields a 128 arcsec x 128 arcsec field-of-view. All ob- 
servations were obtained with a narrow JC'-band filter, Kn, of 
central wavelength 2.165^.m and width 0.33ii-m. The readout 
noise with the double sample readout method used (relative 
sampling) was about 60 e and the dark current < 30 e s~^ 
over most of the array, although there are a number of hot 
pixels. For 5 s integrations (see below) and typical back- 
ground brightness conditions (12.0-13.0 Kn mag arcsec"^), 
our observations are always background limited. 



http: / /www.mso.anu.edu.au/observing/2.3m/CASPIR/ 
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For all observations except biases and standard stars, 
the integration time for a single frame was 60 s, resulting 
from the average of twelve 5 s exposures (individual expo- 
sures were not saved). The total integrations were built-up 
using a script taking 5 frames of each field, each dithered 
by 4" (8 pixels; the biggest bad pixel patch) parallel or per- 
pendicular to the galaxy major-axis, and interspersed by sky 
frames. The lot was bracketed by dark and bias frames and 
repeated typically 7 times, for a total on-source integration 
time of 35 min per field, slightly less on sky. For a galaxy re- 
quiring mosaicing of 2 fields, for example, the sequence was: 
bias, dark, sky, field.l [dith 0], field_2 [dith 0], sky, field.l 
[dith 1], field_2 [dith 1], sky, field.l [dith 2], . . . This ensures 
sufficiently fast sampling of the sky and allows to average 
out bad pixels when the frames are combined. Using the in- 
strument rotator, the CCD rows were always aligned (and 
centered) on the galaxy major-axis, facilitating mosaicing. 
Mosaics of 1 to 5 fields (with 25 per cent overlap) were 
required depending on the galaxy, and all mosaics include 
plenty of sky area around the target. 

The data reduction was mostly carried out within the 
CASPIR IRAF (Image Reduction and Analysis Facility) 
package, following as closely as possible the procedures rec- 
ommended in the CASPIR manual. Dark, sky, and object 
frames were first bias-subtracted and linearized, and the 
skies and objects dark-subtracted. Sky and object frames 
were then flatfielded and the sky subtracted from the ob- 
jects. All the frames of a given object were then combined 
together with proper offsets to form the final image. The 
bias and dark frames used for subtraction were nightly (or 
fractions of a night if deemed necessary) medians of all such 
frames. For sky subtraction, running medians of a few sky 
frames were used, time-centered around the object frame 
being reduced. A fiat was obtained every night from the dif- 
ference of domes frames taken with an incandescent lamp 
on and off. For most objects, a mosaic was created for each 
dither position (e.g. for field_l [dith 0] and field_2 [dith 0] 
in the example above) using blind offsets (the offsets given 
to the telescope), and all such mosaics were then added by 
registering compact sources (galaxy nucleus and foreground 
stars). This procedure minimizes registering errors due to a 
drift of the telescope pointing over time, and is required be- 
cause a guide star common to all fields of a mosaic could not 
be obtained for most targets (the overlap of the 5!2 guider 
fields being minimal or non-existent) . A common background 
across the different fields of each mosaic was imposed by sub- 
tracting the proper constant from each field when necessary 
(imperfect sky subtraction). 

A photometric zero-p oint on the Carter SAAO system 
dCarter fc Meadowj [l995l) was derived nightly by observ- 
ing IRIS (Infrared Imaging Spectrograph) photometric stan- 
dards throughout the night. The i&i-band atmospheric ex- 
tinction correction coefficient was also calculated each night 
by following a few standard stars over a range of airmasses 
(typical value « 0.10 mag airmass"^). The galaxy images 
were also corrected for Galactic extinction, although this is 
negligible for most objects. When frames were obtained un- 
der less than photometric conditions, their fluxes were first 
scaled to the level of photometric ones before being com- 
bined. The final images generally have a background level 
slightly different from zero, and this was again corrected by 
subtracting a constant. 



4 IMAGES 

4.1 Morphological features 

Figures0and|21show, respectively, our data for galaxies with 
a B/PS bulge and for those of the control sample. Each panel 
contains, from top to bottom, first a Digitized Sky Survey 
(DSS) image of the galaxy (for comparison), second our Kn- 
band image and third an unsharp-masked iSji-band image 
(see below), all spatially registered. The surface brightness 
profiles which follow will be discussed in the next section 
(§ lU- As our main interest here is the galaxy bulges, we 
only show the inner parts of our mosaics, which generally 
extend up to and beyond the large-scale discs. 

Comparing the DSS and iih-band images, the advan- 
tages of working in the NIR become obvious. Most impor- 
tantly, dust extinction within the galaxies, when present, 
is drastically reduced (see, e.g., IC 2531 and NGC 7123), 
although it remains non-negligible in a few objects (e.g. 
NGC 1032). Generally, the B/PS and associated features are 
also sharper at _Kn-band, although this is hard to quantify. 
For a similar reason, unsharp-masking, which enhances local 
extrema, is ideal to highlight readily apparent morphological 
features and reveal weaker ones. 

The unsharp-masked images of Figures were ob- 
tained by median-flltering the ifii-band images, that is by 
replacing the value of each pixel by the difference between 
it and that of the median within a centered circular aper- 
ture. The size of the aperture was held fixed across each 
image, but it was varied from galaxy to galaxy to best high- 
light the features of interest here: X shapes, secondary disc 
enhancements, . . . Bigger apertures yield unsharp-masked 
images more similar to the original _KQ-band image, while 
smaller apertures enhance smaller scale features such as cen- 
tral discs, spiral arms and remaining dust lanes. The aper- 
tures used in Figures are thus tightly correlated to the 
characteristic size of the morphological features enhanced. 
We also note that the unsharp-masked images were obtained 
from Ah-band images with the foreground stars interpolated 
over. This process is not perfect, however, and can create 
artefacts in the unsharp-masked images, especially near the 
equatorial plane of the galaxies where the luminosity gradi- 
ents are steep. One must thus be careful when interpreting 
the results (see, e.g., NGC 2788A and NGC 6722). 

As the same morphological features are present in many 
galaxies, we provide a generic description of those features 
below and list the specific features observed in each object in 
columns 2-7 of Table |5| We note however that, for some ob- 
jects, some features are best seen in images unsharp-masked 
on a spatial scale different from that shown in Figures 
To create Table |5| we thus analysed images unsharp-masked 
with circular apertures of radii ranging from 2.5 arcsec (5 
pixels) to 17.5 arcsec (35 pixels). 

OfF-centered X (OX). Many galaxies show an X 
shape structure in the bulge region, usually most easily vis- 
ible in the unsharp-masked ifii-band image. In many cases, 
however, if one follows the branches of the X, the ridges do 
not cross in the centre of the galaxy but rather fall short 
of it. This leads to a '> <' feature rather than a true 'X', 
which we dub an off-centered X (OX). Excellent examples 
are NGC 128 and NGC 6771, but there are many. 

Centered X (CX). In a comparable number of cases, 
the branches of the X do cross in the center of the galaxy. 
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Figure 1. Images and surface brightness profiles of tlie galaxies with a B/PS bulge. From top to bottom, each panel shows first a DSS 
image of the galaxy, second our iiji-band image, third an unsharp-masked iiTn-band image and last major-axis (fainter) and summed 
(brighter) surface brightness profiles, all spatially registered. Contours for the Xn-band image are spaced by 0.5 mag arcsec"^ and the 
faintest countour is indicated in the bottom-right corner of the panel. 
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Figure 1. Continued. 
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Figure 1. Continued. 
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Figure 1. Continued. 



© 0000 RAS, MNRAS 000, 000-000 



K-band observations of boxy bulges. I. 




Position (arcaec) Position (aroaec) 




Figure 1. Continued. 



© 0000 RAS, MNRAS 000, 000-000 



10 M. Bureau et al. 



1 



Iff 

p. 





-» =!*&: "sSO -IB a IS SO 46: .w. 

Foaltloii (aroaee) 




■*m:im; ~m- -X -ib o IB 30 46 . .« 
Positioii (oroaec) 




«m .-rSD -ID 10 ea fjaa' m 
Position (aroaee) 




Poaltlqtn (arcaec). 
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Figure 2. Same as Figure^but for the control sample galaxies. 
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and we dub this feature a centered X (CX). Although the 
transition between off-centered and centered X is somewhat 
ill-defined, good examples of centered X features include 
NGC 1381 and IC 2531. 

Minor-axis extremum (ME). Many galaxies show 
a rather narrow and elongated local maximum along the 
minor-axis, which we dub minor-axis extremum (ME). They 
are generally only visible in 7&i-band images unsharp- 
masked on a smaller scale than that presented in Figuresm^] 
so we show the cases of NGC 1381 and NGC 3203 in FigureOl 
for illustrative purposes. 

Secondary maxima (SM). Many galaxies display lo- 
cal maxima away from the center along the major-axis, which 
we dub secondary maxima (SM). While sometimes readily 
visible in the original i&i-band image (e.g. NGC 4710 and 
ESQ 597-G036), most are easily seen only in the unsharp- 
masked images (e.g. ESQ 151-G004 and NGC 3203). 

Spiral arms (SA). A number of galaxies display nar- 
row and elongated local maxima which are not aligned 
with either axes but rather are typically slightly offset from 
the major-axis, and sometimes even partially wind around 
the bulge forming a pseudo-ring. Those features are gener- 
ally bisymmetric and are again most easily visible in the 
unsharp-masked ifn-band images. We dub them spiral arms 
(SA) for obvious reasons. Good examples include NGC 4469 
and NGC 5746 (but see also the much weaker example of 
NGC 3203 in Fig.EJ. 

In Table 121 we finally give an indication of the inclina- 
tion of the galaxies, labeling them either inclined (I), when 
the galaxy is obviously more than a few degrees from ex- 




PaaltJoa (arcsec) 



Figure 3. Kh-band images of NGC 1381 (top) and NGC 3203 
(bottom) unsharp-masked on a scale of 2.5 arcsec. A long and nar- 
row extremum is clearly seen along the minor-axis of NGC 1381, 
while a weaker one is observed in NGC 3203. We dub both fea- 
tures a minor-axis extremum (ME). NGC 3203 also reveals weak 
bisymmetric elongated maxima slightly offset from the major-axis, 
which we associate with spiral arms (SA) in a nearly but not 
exactly edge-on disc. Those spirals arms are not as apparent in 
images unsharp-masked on larger scales. 



actly edge-on, quasi edge-on (Q), when the galaxy is ap- 
parently away from edge-on by at most a few degrees, and 
exactly edge-on (E), when there is no indication at all that 
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Table 2. Galaxy features. 



Galaxy Images Profiles Notes 

OX CX ME SM SA I MA SU 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



B/PS bulges 



NGC 128 


OX 




ME 


SM 




E 


3 


3 


^Varped disc {in = mode?) 


ESQ 151-G004 


OX 






SM 


SA 


Q 


3F 


3 




NGC 1381 




cx 


ME 






E 


3 


3 


Upwardly curved CX (M)? 


NGC 1596 


OX? 




ME 






E 


2 


2 




NGC 1886 




cx 


ME 


SM 




Q 


3F 


3 




NGC 2310 




cx 




SM 


SA 


I 


3F 


2 


Inner ring with outer spiral arms 


ESO 311-G012 


ox 




ME 


SM 




E 


3 


2 




NGC 2788A 




cx 


ME? 


SM 




Q 


4F 


4 




IC 2531 




cx 




SM 




E 


3F 


2 




NGC 3203 


ox 




ME 


SM 


SA 


Q 


3F 


2 




NGC 3390 




cx 




SM 


SA 


Q 


3F 


2 




NGC 4469 


ox 




ME? 


SM 


SA 


I 


3F 


3 




NGC 4710 


ox 






SM 


SA? 


Q 


4F 


3F 


Inner surface brightness plateau? 


PGC 44931 




cx 


ME? 


SM 


SA 


I 


3F 


2 


Inner ring with outer spiral arms? 


ESO 443-G042 


ox 






SM 




Q 


3F 


3 




NGC 5746 




cx 




SM 


SA 


I 


3F 


3 


Inner ring with outer spiral arms 


IC 4767 


ox 






SM 




Q 


3F 


3 




NGC 6722 


ox 




ME? 


SM 


SA 


I 


3F 


2 


Inner ring with outer spiral arms, warped 


NGC 6771 


ox 




ME 


SM 




Q 


3F 


3 




ESO 185-G053 




cx? 




SM 




E 


3F 


3 


'( ) '-shaped structure 


IC 4937 


ox 






SM 


SA? 


Q 


3F 


2 




ESO 597-G036 


ox 




ME? 


SM 




E 


3F 


3 




IC 5096 




cx 


ME 


SM 




Q 


3 


3 




ESO 240-GOll 




cx? 


ME? 




SA 


I 


3 


2 


Inner ring with outer spiral arms 



Control sample 

NGC 1032 ME E 2 2 

NGC 3957 CX ME SM Q 3F 2 

NGC 4703 ME SM Q 4F 3 

NGC 5084 ME SA? Q 2 2 

NGC 7123 CX ME E 3 2 

IC 5176 ME SA? E 2 2 



Notes: For the images, the features listed were identified on images unsharp-masked on scales ranging from 2.5 to 17.5 arcsec. OX: 
off-centered X feature, CX: centered X feature, ME: minor-axis extremum, SM: secondary maxima along major-axis, SA: spiral arms, I: 
inclination. For the inclination, the galaxies are classified as edge-on (E), quasi-edge-on (Q), or inclined (I). For the surface brightness 
profiles, we list the number of clearly distinct regions. MA: major-axis profile, SU: vertically-summed profile. The letter 'F' denotes a 
fiat intermediate region. 



the galaxy is not exactly edge-on. There is obviously a con- 
tinuum of inclinations, but the first class generally includes 
galaxies where spirals arms are easily visible (e.g. NGC 4469 
and NGC 5746), while the second includes mostly galaxies 
with weak spirals arms or a slightly off-centered (in the ver- 
tical direction) dust lane (e.g. NGC 3203 and NG C 2788A) 
This classification will be especially important in lPaner ll 
dealing with the vertical surface brightness profiles, since 
disc projection effects can not be neglected for discs signifi- 
cantly away from edge-on. 



4.2 Incidence of morphological features 

The most striking result from Figures is that the mor- 
phology of the galaxies with a B/PS bulge is much more 
complex than that of the control sample galaxies, in the 
sense that galaxies with a B/PS bulge contain more of the 
morphological features described in 14.11 (and best-revealed 



by unsharp-masking; see Table Indeed, at least 50 per 
cent of galaxies with a B/PS bulge possess an off-centered 
X feature and at least 38 per cent a centered one, while only 
33 per cent of the control galaxies have either. Similarly, 88 
per cent of the galaxies with a B/PS bulge have secondary 
maxima along the major-axis, while only 33 per cent of the 
control galaxies do. For the presence of spiral arms, the ra- 
tios become, respectively, at least 38 per cent and perhaps 
nil. Only for the minor-axis extrema are those percentages 
inverted, with at most 58 per cent of galaxies with a B/PS 
bulge possessing one but all control galaxies. Those statistics 
thus suggest that both centered and off-centered X features, 
secondary major-axis maxima and spiral arms are prefer- 
entially associated with B/PS bulges, while only minor-axis 
extrema are found preferentially in other bulge types. In fact, 
the contrast between the main and control samples would be 
even greater if, as we will argue in S 16.51 the control sample 
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was not contaminated by weak B/PS bulges (NGC 3957 and 
NGC 4703). 

We also note that, given the large number of secondary 
maxima observed at K-hand (in fact, more than in the opti- 
cal images) , it is unlikely that they are due to the presence of 
obscuring dust filaments inclined with respect to th e equato- 
rial plane, as was suggested by various authors ('e.g. lSandaed 
ll96L:.Whitmore fc BelL.1988.1 . 

4.3 Origin of morphological features and the 
orbital structure of bars 

Although the accretion of exte rnal material can give 
rise to X-shaped featur e s (e.g . iBinnev fc Petroul Il985l : 
iHernauist fc QuinnI Il988l 119891 : IWhitmore fc Belj Il98d) . 
those are generally centered and it is unlikely that any long- 
lasting off-centered X could be produced. Since we observe a 
majority of off-centered X in our sample, and since there is no 
obvious discontinuity between the two, it is more attractive 
to look for a unique mechanism which can simultaneously 
exp lain both centered an d off-centered stable X features. 

lAthanassoulal l|2005ll unsharp-masked a high-quality A^- 
body simulation of a barred disc. The results provide an 
essentially perfect match to the secondary maxima discussed 
above and, depending on the viewing angle to the bar, to 
the off-centered and centered X features, suggesting a link 
between those features and edge-on bars. We discuss below 
the possible origin of this link through the orbital structure 
of barred discs and then discuss some caveats. 

The most important orbit families in 3D bar models are 
those of the xi tree, all elongated parallel to the bar and 
located within corotation. They include the xi family itself 
(restricted to the equatorial plane) and many other families 
bifurcati ng from the xi at vertical res onances of increasing 
energies. ISellwood fc WilkinsonI il993l) provide a general in- 
troduction, but we adopt here the notation of ISkokos et al.l 
|[2002a b) and call them xivi, xiV2, xiv^, etc. The morpho- 
logical features identified in § 14.11 can then be reproduced 
by sup erposing orbit s of the appropriate sha. pes, as done 
bv lPatsis et al.1 ll20o3) . |Patsis fc Xilouri^ tOOdi) recently ex- 
plored very similar ideas, with identical conclusions. 

The centered and off-centered X features, in particular, 
can both arise out of families extending out of the equa- 
torial plane in the vertica l direction. The xiVi fami ly (also 
called banana orbits; e.g. IPfenniger fc Friedlil [l99lh . bifur- 
cating from the xi family at the 2:1 vertical resonance, has 
the largest vertical extent and is shaped like a smile or a 
frown. As shown by the orbit superpositions of lPatsis et al.l 
ll2002i) . xivi orbits seen side-on have a morphology entirely 
consistent with an off-centered X. Depending on the model 
(i.e. mass distribution and pattern speed) and viewing angle, 
centered X features can also be created. 

The X1V3 and X1V4, orbit families (bifurcating from the 
xi family at the 3:1 vertical resonance and shaped like a '~') 
and the X1V5 orbit family (bifurcating at the 4:1 resonance 
and shaped like a 'w') can also give rise to centered or off- 
centered X featu res, depe nding on the model and viewing 
angle (see, again. iPatsis et alJl2003) . This is also the case for 
the zs.is family, which does not bifurcate from the xi orbits 
but is morphologically similar to the X1V4 family, although 
it was present in only one mode l. 

The orbit superpositions of lPatsis et al.l i2002l) further 



show that the orbit families described above (as well as 
higher order families such as X1V7, xivs and xivg) can give 
rise to a number of maxima along the major-axis, similar 
to the secondary maxima identified in § 14.11 Those max- 
ima generally occur at larger radii than the X features and 
near (but within) the ends of the bar, often where orbits 
have loops. An alternative and more straightforward expla- 
nation is that the secondary maxima observed are simply 
the edge-on projections of inner rings (edge brightening), 
known to exist in a large fraction of bar red spiral galax- 
ies and loca ted at the end of the bar (e.g. lKormendvlll979| : 
lButalll995h . Inner rings naturally form in gas-rich discs 
unde r the influence of rotating bars (e.g. Schwarz 1981^, 
Il984l : iBvrd et alJ but th ey also develop in purely 

dissipationless simulations (e.g. ^ Athanassoula fc Misiriotia 
|20b2; Bureau & Athanassoula 2005). Either way, both mech- 
anisms support a relationship to bars, and both mechanisms 
may coexist. Accretion scenarios producing centered X fea- 
tures, on the other hand, will generally not produce sec- 
ondary maxima along the major-axis (even less so system- 
atically located outside of the X feature itsef), making them 
even less attractive. 

The orbit superpositions of IPatsis et alj ll2002h show 
that the barred orbit families considered above can also give 
rise to the minor-axis extrema described in § 14.11 especially 
when seen end-on. The Zi orbits are particul arly interest- 
ing in this respect but were little studied by ISkokos et al.l 
i2002al lbl) and deserve more attention. This is encouraging, 
but the fact that all control sample galaxies also show a 
minor-axis extremum suggests that this feature is not re- 
lated to B/PS bulges (and thus presumably to bars). We do 
not properly understand the origin of this feature. 

As we mentioned in § 14.11 because of the morphological 
similarities, we interpret the local maxima elongated parallel 
to but slightly offset from the major-axis of the galaxies as 
spiral arms. If this association is true, which we have no rea- 
son to doubt, then obviously spiral arms can only be detected 
in galaxies that are not perfectly edge-on. In fact, we have 
used the presence of spirals to quantify the inclination of 
the galaxies in Table |5| Since we have no other independent 
and reliable way to estimate the inclination of the galaxies 
precisely (as we do not know the intrinsic thickness of the 
discs), we can not guarantee that the inclination distribu- 
tion is similar between the sample of galaxies with a B/PS 
bulge and the control sample. Thus, a reliable comparison of 
the incidence of spiral arms between the two samples is not 
possible, and we do not discuss this feature at length. 

Suffice it to say that the apparently slightly higher inci- 
dence of spiral arms in galaxies with a B/PS bulge is consis- 
tent with the suggestion that B/PS bulges are simply thick 
bars viewed edge-on. Indeed, bars are very good at driv- 
ing (grand design) spiral patterns (e.g. j Athanassoula 1980; 
iThielheim fc Wo lff 1984: Elmegreen fc E lincgrccn 1989) and 
spiral arms are preferentially found in barred ga laxies (e.g. 
lElmegreen fc ElmegreenI Il989l: HutaeTaO Eo^ . Further- 
more, for the few galaxies where a (partial) inner ring with 
outer spiral arms is visible (e.g. NGC 5746), the ring is sys- 
tematically located just at the end of the flat intermediate re- 
gion of the major-axis surface brightness profiles (see §|3 all 
such galaxies have an intermediate region, and it is flat in all 
cases but one). If, as we argue in §|^ the end of the intermedi- 
ate surface brightness profile region also marks the end of the 
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bar, then this trend is entirely consistent with observations 
of mo re facc-on barred spiral galaxies (e.g. lKorme ndvlll979|; 
&i^y-995) an d the expectat ions from theory (e.g. Schwar a 
1981 n 9841: IRvrd et a,lJll 9'9l. 

On a different note, we must point out that while the 
(barred) picture proposed above is attractive and supported 
by much independent evidence (see § 0, it may not be 
unique. Indeed, as simple tests show, surface brightness dis- 
tributions created from nested rectangular isophotes can 
also give rise to centered and off-centered X-shapes when 
unsharp-masked, independently of their origin. As axisym- 
metric box y distributions can in principle also exist (e^. 
Binnev&Petrou, .1985: May, van Albada fc Norm an 19851 : 
RowlevlllQsir it is possible that bars are not uniquely re- 
lated to the complex features observed. However, as the N- 
body simulations and orbit properties discussed above show, 
barred models are consistent with the unsharp-masked fea- 
tures in much details, while this remains to be shown for 
other models. 



5 SURFACE BRIGHTNESS PROFILES 
5.1 Profile features 

In addition to the DSS, _Kh-band and unsharp-masked Kn- 
band images, the panels of Figures also contain two sur- 
face brightness profiles. The fainter ones represent the major- 
axis surface brightness profiles, extracted using a constant 
position angle, an approximation that is valid in the inner 
parts of all galaxies and in the outer parts of all but a few 
(e.g. NGC 128, NGC 6722, NGC 6771 and ESQ 597-G036). 
The brighter surface brightness profiles were obtained by 
summing the data in the vertical direction, as if the galax- 
ies were infinitely thin, again assuming a constant position 
angle for the major-axis. In both cases, the /(h-band images 
with the foreground sources interpolated over were used, and 
the vertical profiles were summed until the noise level of the 
image was reached. While the noise level varies slightly from 
object to object, such a limit is sufficient for our purposes. 

While our summed surface brightness profiles are not 
strictly equivalent to azimuthally-averaged surface bright- 
ness profiles in more face-on galaxies, or to cuts at a given 
position angle (i.e. they are also integrated along the line- 
of-sight through the edge-on discs), a comparison is never- 
theless possible and useful. Indeed, because the luminosity 
of discs generally decreases rapidly with radius, the light at 
any given position in our profiles is still dominated by that 
near the same (cylindrical) r adius. O u r profi les can also be 
compared to those in Liitticke et al who showed 

profiles (cuts) parallel to but offset from the major-axis for 
a number of edge-on spiral galaxies. 

In columns 8 and 9, Table |5| lists the number of dis- 
tinct regions present in the surface brightness profiles of each 
galaxy, each region being typically separated by a clear ra- 
dial break (although one must be cautious with badly in- 
terpolated foreground sources). From the surface brightness 
profiles of face-on galaxies, we would generally expect each 
profile to show only two distinct regions if the galaxy is ax- 
isymmetric. A first region at small radii, generally associated 
with the bulge, where the profile is steep (generally steeper 
than an exponential profile), and a second region at larger 
radii, generally associated with the disc, where the profile is 
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Figure 4. Cartoon description of our surface brightness profiles, 
as tabulated in Table|^ Two, three, and four-region surface bright- 
ness profiles are shown, with or without a flat intermediate region 
(denoted by the letter 'F'). Note that no '2F' profile was actually 
observed. 

exponential^. However, such two-region profiles are rare in 
our data, especially along the major-axis (but see NGC 1596 
and NGC 1032). Most profiles show an additional region at 
intermediate radii, where the profile is shallower than the 
outer disc, even often flat or slightly rising with radius (de- 
noted by the letter 'F' in Tabled e.g. ESO 151-G004 and 
NGC 2310). Such three-region surface brightness profiles are 
often referred t o as Freeman Type II profiles in the literature 
(lFreemaJl97d) . A few galaxies also show two distinct inter- 
mediate regions, one shallow following the steep central com- 
ponent and one steep preceding the outer exponential disc, 
for a total of four regions (e.g. NGC 2788A and NGC4710). 
We do not discuss (or count) here apparent radial breaks in 
the outer discs, since they are presumably only weakly re- 
lated to the bulge structure. Figure|l|shows a cartoon version 
of the various types of profiles identified in Table |21 

5.2 Incidence of profile features 

Again, the most obvious and basic result from Figures 012] 
is that the surface brightness profiles of the galaxies with a 
B/PS bulge are much more complex than those of the con- 
trol sample galaxies, especially along the major-axis, in the 
sense that the surface brightness profiles of galaxies with 

^ We remind the reader that the surface brightness profile of an 
edge-on exponential disc is not exactly linear in a (logarithmic) 
surface brightness vs. projected radius plot, but slightly curved 
(see, e.g.. |Pohlen et al...200a) . 
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a B/PS bulge typically contain a larger number of distinct 
regions separated by clear radial breaks (see Table I^J . In- 
deed, along the major-axis, 96 per cent of the galaxies with 
a B/PS bulge have a surface brightness profile with 3 or more 
regions (Freeman Type II profile), 78 per cent of which have 
a flat intermediate region, while the corresponding fractions 
for the control sample are, respectively, only 50 and 66 per 
cent. The statistics are not as sharply contrasted when the 
profiles summed along the vertical axis are considered, but 
the differences remain significant. For example, 58 per cent 
of the galaxies with a B/PS bulge have a summed surface 
brightness profile with a least 3 regions, while only 17 per 
cent of the control gala xies do. Those fractions are in rough 
agreement with those of lLiitticke et all i2000lj) , who carried 
out a similar analysis. As for the morphological features, 
it thus appears that complex surface brightness profiles are 
preferentially associated with B/PS bulges. And again, the 
contrast between the main and control samples would be 
sharper if the control sample had been better selected (i.e. 
if NGC 3957 and NGC 4703 has been excluded; see Sl^. 

5.3 Origin of profile features and relationsiiip to 
bars 

The fact that practically all galaxies with a B/PS bulge have 
a major-axis surface brightness profile with an intermediate, 
often flat region is particularly important, as it is inconsis- 
tent with an axisymmetric bulge and exponential disc model, 
no matter how steep the bulge proflle is (i.e. no matter what 
the Sersic n index of the bulge is). There is thus apparently 
a third photometric or morphological 'component' dominat- 
ing at intermediate radii. We argue below, however, that in 
such Freeman Type II proflles, both the steep inner region 
and the flat intermediate region are ultimately caused by a 
single bar viewed edge-on, without the necessity of a classic 

bulge; 

iBureau fc Athanassoulal JioO^ followed the temporal 
evolution of the major-axis surface brightness profile in a 
large number of barred A''-body simulations viewed edge- 
on. They showed that bar formation and evolution within 
an (initially) exponential disc is associated with the buildup 
and continued growth of a dense central region, resulting 
in a central peak (which would normally be identified with 
a bulge) and with the formation and gradual flattening of 
an intermedi ate region, in addition to the outer exponential 
disc (see also lAthanassoula fc Misirioti3l2002t lAthanassoulal 
I2OO3I) . The central peak extends to 1-1.5 original disc scale- 
lengths while the intermediate region extends to the end of 
the bar, well beyond the central peak and the thickest part 
of the bar. The existence of three surface brightness pro- 
file regions is thus totally normal (and expected) in barred 
galaxies, but it remains unexplained in classic bulge models, 
unless there is an additional bar. In fact, four-region surface 
brightness pr ofiles are also often produced in barred models 
(e.g. Fig. 3 in lBureau fc Athanassoulall2005l) . 

As expected from the elongated boxy/peanut shape of 
the bar, the ratio of the length of the thickest part of the bar 
(or that of the central peak) to that of the flat intermediate 
region increases in A^-body simulations as the viewing angle 
increases (from end-on to side-on) . This trend is also present 
in our data, where the ratio of the length of the thickest part 
of the bulges (or that of the central peak) to that of the in- 



termediate region is generally larger in peanut-shaped bulges 
than in boxy ones (compare, e.g., NGC 5746 and NGC 6722 
to NGC 1381 and IC 5096). There is however much variety, 
as many strong peanut-shaped bulges have a small ratio (see, 
e.g., NGC 1886 and ESQ 443-G042), and likely many causes 
for it. First, as far as the surface brightness profile and bar 
shape are concerned, there is a certain degeneracy between 
bar strength and viewing angle (a strong bar seen at interme- 
diate viewing angle may appear similar to a weake r bar seen 
exactly side-on; e.g. iBureau fc Athanassoulall2005ll . Second, 
the ratio described evolves in A^-body simulations, gener- 
ally decreasing with time as the bar strengthens and length- 
ens but the central compone nt changes little (see, again, 
IBureau fc AthanassoulallioOSl'l . Third, gas inflow and subse- 
quent star formation (absent from the pure A-body models) 
may substan tially modify the light distribution in the centra l 
regions fe.g. iHeller fc Shlosmai]ll994l:lFriedh fc Benzlll995l) . 
Fourth, the presence (or absence) of rings and lenses will 
co mplicate the comparison . We a lso note that as the profiles 
of IBureau fc Athanassoulal ||200^ were integrated vertically 
over 0.5 original disc scalelengths, they are somewhat inter- 
mediate between our major-axis and summed proflles. 

In fact, the variety observed in the ratios of the length 
of the thickest part of the bulge (or that of the central peak) 
to that of the intermediate region may well be dominated by 
the range of bar strengths present in the sample, rather than 
by the range of viewing angles . This is also suggested by the 
data of lLiitticke et all i2000lJl . and it would similarly imply 
that the observed boxy-peanut 'sequence' is not mainly one 
of viewing angle, but is instead driven by a varying (i.e. 
increasing) bar strength. The small ratios observed in strong 
peanut-shaped bulges such as NGC 1886 and ESQ 443-G042 
can then be explained only if the central surface brightness 
peak and the thick part of the bulges are shorter in stronger 
bars (in a relative sense). This is natural if the inner part 
of the surface brightness proflles is dominat ed by a disc- 
like bu lge, as will be argued in § 16.41 Indeed, lAthanassoulal 
(Il992al lbl) showed that the size of (bar-driven) nuclear gas 
discs is limited by the outer inner Lindblad resonance (i.e. by 
the outermost X2 orbit, elongated perpendicular to the bar), 
which itself inversely depends on the bar strength. Then, the 
stronger the bar, the smaller the nuclear disc. 



6 DISCUSSION 

6.1 Radial redistribution of material 

lAthanassoulal (|2002| . |2003f) argued that much of the bar- 
driven evolution in discs is due to a transfer of angular mo- 
mentum from the inner (barred) disc to the outer disc (be- 
yond corotation) and the dark halo, leading to a strength- 
ening, lengthening and slowing down of the bar with time. 
To isolate the effects of angular momentum exchange, which 
leads to a radial redistribution of matter, we must preferen- 
tially consider the vertically-summed surface brightnes pro- 
files of Figures and |5| While fewer galaxies with a B/PS 
bulge exhibit 3 or more regions in the summed proflle than 
along the major-axis, a majority still do, while only one of 
the control sample galaxies does (see Table |5J. 

In scenarios where the formation and evolution of discs 
is dominated either by collapse or the gradual accretion 
of discrete components (e.g. cannibalized dwarf or satellite 
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galcixies), there is no obvious mechanism to create the radial 
breaks in the disc density distribution observed in our sam- 
ple galaxies, even less to explain the spatial correlation of 
those breaks with the ionized-gas and stellar kine matics (see 
iBureau fc Freemanlll999l:IChung fc BureaulEooi) . Similarly, 
the traditional interpretation of Freeman Type II profiles as 
inwardly truncated discs neither proposes an origin for the 
truncation nor explains why the truncation is only partial 
and the surface brightness profiles systematically have the 
specific shape observed. In bar-driven evolution scenarios, 
however, a break naturally occurs at the end of the bar (e.g . 
lAthanassoula & Misiriotis" 200^ lAthanassouTal l2002l 12003: 
iBjireau fc Athanassoula 2005). If this is indeed the case here, 
then the break in a surface brightness profiles at the end of 
the intermediate region should systematically mark the bar's 
end. 

Our major-ax is surface brightness p rofiles match 
those of A tliMiassoiilafc Misirioti'3 (|20o3) and 
iBureau fc Athanassoula |2005r ! and as we mentioned 



already the kinematics of the sample 
well with barred galaxy models (see 


galaxies agree very 


IBureau fc FrccmarJ 


I2OO4: 


'Athanassoula fc Bureau' 'l99d: 


IChune fc Burea^l 




.Bureau fc Athanassoula .2005). A Quick comparison 



further shows that the position where the rotation curve 
flattens, ge nerally associated with the end of the bar (e.g. 
IBureau fc Athanassoula 2005), always occurs very near 
(although generally slightly within) the break in the Kn- 
band surface brightness profile. Thus, our surface brightness 
profiles are consistent with those expected from bar-driven 
evolution scenarios, and the complex surface brightness 
profiles observed appear to be a direct consequence of the 
transfer of angular momentum and the radial (and vertical) 
rearrangement of material by the bar. 

Most barred galaxies appear to have fast bars, 
whereby corotation, defined as the radius where the stars 
and bar pattern rotate at the same speed (and mark- 
ing the transition between material loosing and gain- 
ing angular momentum), is located just beyond the 
end of the bar (see, e.g.. iGers sen. Ku iiken fc Merrifieldl 
I2OO3I: lAguerri. Debattista fc Corsinil l2003l) . Since the pat- 
tern speeds of the bars in our sample are unknown, a di- 
rect test of this is impossible. However, for the few galax- 
ies with an apparent inner ring (see Table (21 , the ring ra- 
dius is systematically (roughly) equal to that of the break at 
the end of the intermediate surface brightness region. Given 
that inner rings are generally thought to occur near the 
inner 4:1 (ultra-harmonic) and corotation resonances (e.g . 
Schwarzlll98ll. I l984': 'Athanassoula et ahl ligsl lButalll995l: 
Patsis. Skokos fc Athanassoula 2003b), our galaxies are con- 



sistent with harbouring fast bars. 

It is also interesting to note that the two latest galax- 
ies in our sample, ESO 240-GOll and IC 5176 (with very 
small bulges and thus presumably very weak bars, if any) 
both have completely featureless (single) exponential outer 
discs. They thus offer useful benchmarks against which more 
complex systems can be compared. 



6.2 Vertical redistribution of material 

To isolate the effects of the vertical redistribution of material 
within bars, as predicted by models where the bar buckles 
and thickens (e.g. 'Combos fc Sa nders, .198L : .Combes et al.i 



ll990l : lRaha et al.ll99 j) . we compare here the major-axis and 
vertically-summed surface brightness profiles of our sample 
galaxies. 

A cursory examination of our major-axis and summed 
surface brightness profiles shows that, for most galaxies with 
a B/PS bulge, they differ markedly from one another. As 
already pointed out, the vertically-summed profiles generally 
show fewer of the features described in §|^and they approach 
more closely a simple bulge plus exponential disc description, 
while the major-axis profiles generally show a steeper inner 
peak and a fiatter intermediate region (the outer regions have 
a comparable slope in both profiles). Excellent examples of 
this are ES0151-G 004 and PGC 44931, but there are many. 

If the scaleheight of the stars in a galaxy is constant with 
radius, one expects the major-axis and vertically-summed 
profiles to have the same functional form but different zero- 
points, i.e. be offset but parallel to one another (although 
not exactly, since the profiles are vertically-summed to a 
given surface brightness level rather than to a fixed num- 
ber of scaleheights) . This is not the case for most galaxies, 
however, even when the major-axis surface brightness pro- 
file is fiat (when summing to a given surface brightness level 
is equivalent to summing to a fixed number of scaleheights) , 
suggesting that the stellar scaleheights are not constant. See, 
again, ESO 151-G004 and PGC 44931, ahhough any galaxy 
with a fiat intermediate profile will do. 

By discussing a single scaleheight, even in regions where 
the inner peak of the surface brightness profiles is non- 
negligible, we have however implicitely amalgamated the 
bulge and disc (no matter how the former is defined) . In the 
classic bulge and disc model, where the bulge and disc rep- 
resent structurally and kinematically distinct compoments, 
this will necessarily lead to a variation of the scaleheight, al- 
though one would naively expect this variation to be mono- 
tonic. The major-axis and vertically-summed surface bright- 
ness profiles of Figure Q however, show that the variations 
are not monotonic. Furthermore, the functional difference 
between the two profiles is systematically greatest in the fiat 
intermediate region, which would never be associated with 
a bulge (even in classic models) and is clearly dominated by 
disc material. Our major-axis and vertically-summed sur- 
face brightness profiles thus show that the radial scaleheight 
variations are real and that they occur in the discs (i.e. they 
are not due to bulge material), in direct contradiction to 
the long accepted wisdom that disc scaleheights ar e radia ll;y 



tne long acceptea wisaom tnat gisc scaieneignts ar e raaia iiy 
constant (see, e.g. . Iron der Kruit fc Searlelll981alin.ll982alB 



IShaw fc Gilmorelll99d : Ide Griis fc van der Kruitlll99d) . In 

fact, our results support a picture where the bulge and disc 
are not intrinsically distinct (either structurally or kinemat- 
ically), but both instead emerge from the rapid radial varia- 
tion of the scaleheight of the disc material (due to bar-related 
vertical r esonances). This will be discussed in more depth in 
IPaper l3. where it will also be shown that the variations are 
as expected from barred Ai'-body models. Viewed simplisti- 
cally, our observations argue for the reassignment of (most 
of) the bulge material to the disc (a B/PS bulge), but this 
is an important conceptual change. 

To our knowledge, a variable stellar scaleheight has only 
been identified and discussed in a few galax ies. In the outer 
parts of NGC 3115. ICapacciolT et al.l (|l98i) assign the vari- 
ation to the end of the disc self-gra vity, while in the inner 
parts of the Milky Wav lKent. Dame fc Fazig (>1991il assign 
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it to a variable ratio of young to old stars (their data in fact 
already show strong hints of the peanut shape of the Galac- 
tic bu lge and a Freeman Type II profile) . Ide Griis fc Pele tieij 
do note an increase of the scaleheight with (projected) 
radius in a large number of edge-on spiral galaxies, but they 
assign the increase to the presence of a thick disc with scale- 
length and scaleheight larger than those of the dominant 
(thin) disc component. Upon closer inspection, a number of 
other studies arguing for a constant scaleheight as a func- 
tion of (projected) radius actually do detect variations, but 
those are generally monotonic and are argued to be consis- 
tent with the infiuence of a structurally distinct bulge (in 
the center) or a thickening of the outer disc (near the edge 
of the optical disc). The scaleheight variations are thus re- 
jected as artefacts unrelated to the domin ant disc, presumed 
thin (e.g. Ide Griis fc van der KruidlToQ^ . While this is de- 
fendable in some cases, it is doubtful in others, and the non- 
monotonic variations observed here and in .Paper III in the 
inner and intermediate regions of our sample galaxies can 
not be explained away in such manners. Possible variations 
in the outer parts of our galaxies (e.g. fiaring) are not of 
direct interest to the bulges' structure. 

6.3 Classic versus 'pseudo'-bulges 

As discussed also in lAthanassoulal ^|2005^ ■ galaxy bulges have 
traditionally been defined either as 1) the steep central com- 
ponent of the surface brightness profile (generally steeper 
and brighter than the inward extrapolation of the o uter ex- 
ponential disc; e.g. ICarollo. Ferguson fc WvseHlQQgT) . 2) the 
thick galactic component (in isophotal terms, clearly stick- 
ing out of the equatorial plane in nearly edge-on objects) 
or 3) the kinematically hot component (the central peak in 
the velocity dispersion profile) . Those three definitions have 
been used interchangeably and have generally been consid- 
ered equivalent and consistent. But is this really the case? 

Many diverse resu lts show those views to be g rossly 
oversimplified (see, e.g.. lWv~ Gilmore fc Franxll997l for an 
early review). For example, there is now much evidence that 
some bulges are really concentrated (thin) discs, or equiva- 
lently that the central parts of discs can have a density pro- 
file st eepe r than that of the outer disc (s ee, e.g.. iKormendvl 
Il993l and iKormendv fc Kennicutl |2004| for reviews). The 
central surface brightness profiles o f most bulges also rarely 
approach the classi c i?^/* law (e.g. lAndredakis et alj[T995l 
iBalcells et ai]l2003l) . Although our own surface brightness 
profiles are not ideal to tackle this last issue, because of the 
line-of-sight integration through the outer disc, many ob- 
jects are nearly exponential in t he center and we specifically 
study the minor-axis profiles in iPaoer ll Most importantly, 
our sample is ideal to test the equivalence of thick structures 
with steep light profiles and high velocity dispersions. 

The images and surface brightness profiles of Figure 
reveal that, in some galaxies, the thick central structure (de- 
fined isophotally) extends past the steep inner region of the 
surface brightness profile (both along the major-axis and ver- 
tically summed). Those two definitions of a bulge are thus 
neither equivalent nor consistent, at least in those cases. 
Our observations thus show that the steep inner region of 
the surface brightness profile is often contained within the 
thick central component. This was already hinted at by the 
poorer optical profiles of iChung fc Bureau. (■2004i'l . although 



iLiitticke et al ] l|2000b) found them to be generally equal. Ex- 
amples include ESQ 151-G004, NGC 2788A and ESQ 443- 
G042, but there are many. In fact, our observations suggest 
that those two bulge definitions only seem to agree systemat- 
ically in galaxies with a boxy-shaped (or round) bulge, such 
as NGC 1381 and IC 5096, while the disagreement is worst 
in some of the the strongest peanut-shape bulges. There is 
no straightforward explanation for this in classic bulge for- 
mation scenarios, but it is a natural consequence of bar- 
buckling/thickening mechanisms, where it is related to the 
bar viewing angle (seen respectively end-on and side-on). 

Having argued that the steep inner region of the surface 
brightness profiles is generally contained within the thick 
part of the galaxies, it is also true that the intermediate 
and often flat region of the surface brightness profiles al- 
ways extends beyond the thick component. Some of the best 
examples are found in the galaxies listed above, but again 
there are many (see, e.g., NGC 1886 and NGC 4469). This 
behaviour is hard to reconcile with classic bulge formation 
scenarios, but it arises naturally in bar-buckling/thickening 
ones. Indeed, while the (projected) surface brightness profile 
at intermediate radii remains shallow or flat within most of 
the barred region of a galaxy (i.e. over the entire depro- 
jected bar length, no matter what the viewing angle is), 
only part of the bar is thick. So, while the fact that most 
thick components are shorter than the flat intermediate re- 
gion of the surface brightness profile is simply a consequence 
that most bars are not seen exactly side-on, the fact that 
all thick components are shorter than the flat intermedi- 
ate region reflects a fact much less widely appreciated, i.e. 
that only a fraction of the bar is actually thick. This point 
has not been stressed in the literature , but it is nice l y illus - 
trated and discussed in some depth in I Athanassoulal (1200511 . 
It is explained by t he orbital structure o f barred discs (e.g. 
ISkokos et aT]l2002allhl: IPatsis et al1l2002ll and is in fact vis- 
ible in most published A^-body models fe.g. iCombe s et ahl 
Il990l : lAthanassoula fc Misiriotyl2002l : lAthanassoulall2005^ . 

6.4 Three-dimensional structure of B/PS bulges 

A comparison of the major-axis and vertically-summed sur- 
face brightness profiles reveals another interesting fact, i.e. 
that the steep inner region of the surface brightness profiles is 
systematically more pronounced and domina nt in the major- 
axis pr ofiles than in the summed ones (see also lLiitticke et ciLI 
l2000bl) . This is a key indication of the 3D structure of galax- 
ies with a B/PS bulge. It shows that most of the material 
found at high z belongs to the flat intermediate region of 
the surface brightness profiles rather than to the inner steep 
component, and that the latter has a smaller characteris- 
tic height than the former (this will be shown more trans- 
parently in Paocr II). Thus, the steep inner region of the 
surface brightness profiles truly appears to be a thin concen- 
trated disc, while the shallow intermediate region appears to 
be thick (and, according to the aforementioned arguments, 
barred). Figure |^ presents a cartoon version of this model, 
con trasted to the classic one. 

lAthanassoulal i2005l) recently proposed a new nomencla- 
ture to differentiate different types of bulges based on their 
formation mechanism. In addition to classic bulges, which 
presumably form rapidly through either collapse or merging, 
she identifies B/PS bulges, which are simply (part of) thick 
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Figure 5. Cartoon description of our B/PS bulge model (left) and the classic bulge model (right). From top to bottom, each panel shows 
the face-on view of the model, the edge-on (side-on) view and the edge-on surface brightness profile, all spatially registered. Galaxies 
with a classic bulge are generally viewed as having both a thick and steep inner (classic bulge) component within a thin exponential disc. 
We suggest that galaxies with a B/PS bulge instead harbour a thin and steep inner disc component (disc-like bulge) within a thick and 
shallow bar (B/PS bulge), followed by an outer exponential disc. 



bars, and disc-like bulges, formed by (bar-driven) gaseous in- 
flow and subsequent star formation. As more than one type 
of bulge can coexist in a single object, the latter two types 
are well suited to a description of the observations presented 
here. 

As noted abo v e. [ Kormendvl (Il993l) and 
iKormendv fc KennicuttI (|2004|) review much of the ev- 
idence for the interpretation of some bulges as concentrated 
discs, which they include in a broader category of 'pseudo'- 
bulges, generally formed secularly (i.e. slowly). This vision 
of disc-like bulges does not explicitely include the extended 
bulges discussed in this paper, nor the many bulges which 
are isophotally thick and kinematically hot. The thick bar 



model advocated in lAthanassoulal ll2005h and here ad- 
dressess all three points simultaneously. In this respect, we 
know that the radial component of the velocity dispersion 
of the Milky Way bu lge is similar to that of its inner disc 
(» 100 km s~^; e.g. [ Soaenhauer. Jones fc Whitford 1992i ; 
iLewis fc Freemanlll989^ ■ and that the velocitv dispersion of 
the galactic bar/bulge is nearly isotropic. So although little 
is known about either, a difference between the vertical ve- 
locity dispersion of the inner disc and that of the bar/bulge 
must underlie their large structural differences. While 
it is unlikely that concentrate d discs are hot vertically , 
(thick) bars appear to be (e.g. iGadotti fc de Souzal I2OO5I : 
■Debattista ot aL.2005.) . 
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One could argue that our model comprises a number 
of distinct building blocks, like the classic models, but those 
blocks are very different and tightly intertwined dynamically. 
First, the large-s cale bar lea ds to th e formation of a concen- 
trated disc, as in|Ko^Mn^3s lll993ll picture, but this disc is 
thin, largely decoupled from the large-scale bar, and only ad- 
dresses the first of the bulge definitions (steep surface bright- 
ness profile). Second, the large-scale bar itself is thick over 
most but not all its length, and it addresses the second defi- 
nition of bulges (isophotally thick structure) . As for the third 
definition of bulges (kinematically hot), it is unclear which 
of the inner disc or large-scale bar dominates the large veloc- 
ity dispersions observed by fChung fc Bu reau ( 2004) in the 
equatorial plane, b ut the models suggest that th e culprit is 
slightly elongated iBureau fc Athanassoulalf2005l) . As men- 
tioned above, face-on ba rs do appear to have large vertical 
velocity dispersions (e.g. iGadotti fc de Souzall2005h . 

The 3D struc ture proposed for B/PS bulges in 
lAthanassoulal (|200^ and here is thus diametrically opposed 
to the classic one. Instead of a thick and steep inner (classic 
bulge) component within a thin exponential disc, we propose 
a thin and steep inner disc component (a disc-like bulge) 
within a thick and shallow bar (a B/PS bulge), followed by 
an outer exponential disc (see Fig. [SJ. Furthermore, while 
the two components of classic bulges are generally thought 
to have formed separately through vastly difi'erent mecha- 
nisms and to interact little with each other, the three 'com- 
ponents' of galaxies with B/PS bulges are largely made of 
the same material, which has evolved and been shaped by 
the single dominant infiuence of the bar, through the weak 
but relentless action of both radial and vertical resonances. 
We do not appear to need additional classic bulges to ex- 
plain our observations, although it is difficult to rule them 
out completely at the very center of our galaxies (i.e. on a 
scale much smaller than the isophotally thick part). 

We note that while this may seem to imply that no stel- 
lar population gradient should exist in barred galaxies, this is 
not true. The radial resonances, once established, will clearly 
affect the behaviour of the gaseous (dissipative) material, 
which in turns feeds star formation. It is thus expected that 
a certain segregration of the stellar populations (mainly due 
to influence of the radial resonances) will develop over time. 
While pure A/'-body models do not capture this richness, 
others models including gas dynamics and star formation 
do (e.g.lFViedh. Benz fc Kennicuttll l994': 'Heller fc Shlosman 
| 1994|: iFriedli fc Bendll995l : but see also iNakasato fc Nomotoi 
I2OO3). However, the current simplicity of those models and 
the ignorance of much of the physics involved (primarily 
regarding star formation) prevent a detailed comparison 
with cu rrent obsvervations, which are in any case limited 
(but seelMartin fc Rovlll994lZaritskv. Kennicutt fc Huchral 
11994: iFisher. Franx fc Illingworthlll99d) . There is no doubt, 
however, that the next step for bar-driven bulge for- 
mation models is to include realistic star formation 
and stellar population information ( chemodynamics; see , 
e.g., 'Samland. Hensler fc Thei^ Il997l : ISpurzem et"ai] 120041: 
iMichcl-Dansac fc Wozniakl booil . 

The fact that the maximal radial extent of the inner 
disc appears to be systematically equal to or slightly smaller 
than the radius where the bar is thickest (see Fig. may 
reflect some form of 'self-regulating' dynamics, i.e. that the 
inner disc is limited by a radial resonance (perhaps the radial 



inner Lindblad resonance, ILR) coinciding with a vertical 
resonance (p erhaps the vertical ILR), as advocated e.g. by 
ICombes et aL (.1990.') . This should be explored further. 

6.5 Control sample 

Due to its limited size compared to the number of galaxies 
with a B/PS bulge in our sample, the control sample was 
always meant more for qualitative checks than for statisti- 
cally robust comparisons. However, it is now clear that some 
galaxie s do not belong in it. For example, Liitticko et aj] 
]|oro3) classified the bulge of NGC 3957 as close to boxy- 
shape d, iBureau fc FreemanI il999l) and IChung fc BureaiJ 
||2003) both found kinematic indications of weak ba rs in 
NGC 3957 and NGC 4703, and lChung fc Bureaul (j^Ooi) also 
found indication of non-axisymmetric motions in NGC 7123. 
The current K-hand imaging also shows centered X fea- 
tures in NGC 3957 and NGC 7123 and secondary maxima 
in NGC 3957 and NGC 4703. AU three galaxies show at 
least 3 regions in their surface brightness profile. The con- 
trast between the B/PS bulge and control samples would 
thus be greater if NGC 3957 and NGC 4703 (and perhaps 
NGC 7123) were appropriately classified as B/PS. 

Even if bars are peanut-shaped, a fraction of round 
bulges should still be expected to be bars, as the latter ap- 
pear round when seen end-on. This fraction is however ex- 
pected to be small, since already for an angle of 10-30° be- 
tween the bar major-axis and the line-of-sight a boxy shape 
is observed (e.g. C ombes ot al. 1990; Liitticko et al. 20Q0b. : 
lAthanassoulal lioOSf) . The large contamination of our con- 
trol sample thus probably arises from the way it was se- 
lected. Indeed, no catalog of edge-on galaxies with a round 
bulge existed when it was constructed, so the control sam- 
ple galaxies were selected from apparently misclassifi ed ob- 
jects in catalogs of galaxies with a B/PS bulge jjarvisl 
119861 : Ide Souza fc dos Anio3ll987l:IShawlll987l) and from the 
Flat Galaxv Catalog "oF lKara^hentse^e^i^^ Carry- 
ing a study along the lines of the current on e and those of 
iBureau fc Freeman (199a) and lChung fc Burea u (2004), but 
for an enlarged and properly selected sample of round bulges 
(for ex ample using the more recent catalog of lLiitticke et al.) 
|2000a|) ■ would thus be valuable. While this represents a sig- 
nificant amount of work, showing with more certainty that 
round bulges behave differently from B/PS bulges would 
greatly strengthen our conclusions. 



7 SUMMARY AND CONCLUSIONS 

We have presented _Kh-band imaging observations of a sam- 
ple of 30 edge-on spiral galaxies, most of which harbour a 
boxy or peanut-shaped (B/PS) bulge. Those data are min- 
imally affected by dust and best trace population II stars, 
where most of the luminous mass resides. Our multi-faceted 
analysis suggests that B/PS bulges are simply the thick part 
of bars viewed edge-on (see Figures and I^J. 

Galaxies with a B/PS bulge tend to have a more com- 
plex morphology than galaxies with other bulge types, more 
often showing centered or off-centered X structures, sec- 
ondary maxima along the major-axis, and spiral-like struc- 
tures. Best revealed by unsharp-masking, those features are 
also observed in three-dimensional N-body simulations of 
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barred discs ijAthanassoulalli'oO.'Th . and can be explained b y 
the orbital structure of bars (see, e.g., IPatsis et alJl20o3) . 
although they need not be uniquely related to them. Only 
minor-axis extrema may be preferentially related to other 
bulge types. Whether taken along the major-axis or summed 
vertically (to simulate a flat galaxy), the surface brightness 
profiles of galaxies with a B/PS bulge are also more com- 
plex, more often showing 3 or more clearly separated regions, 
including a rather shallow or flat intermediate region (see 
Figure 1^. Such Fre eman Type II profiles are expe cted from 
barred galaxies (e.g. lBureau fc Athanassoulal2005l) . but they 
do not have a natural self-consistent explanation in classic 
bulge formation scenarios. 

The radial breaks observed in the vertically-summed 
profiles of our objects provide further evidence of the transfer 
of angular momentum and radial redistributio n of disc mate- 
rial mediate d by the (presumed) bars (see, e.g.. lAthanassouial 
I2OO2I . l2003l) . Furthermore, the spatial correlations of the 
radial breaks with the ionized-gas and stellar kinematics 
|Burcau & Freeman 1999; Chung & Bureau 2004) arc as ex- 
pected for fast bars, currently favoured by observations. The 
differences between the major-axis and vertically-summed 
profiles provide evidence for abrupt variations of the scale- 
height of the disc material. This is, again, as expected from 
the diverse orbital families and vertical resonances and insta- 
bilities present in barred discs, but contrary to conventional 
wisdom. A quantitative and robust analysis of those scale- 
height variations and a comparison with A''-body simulations 
will appear in future papers of this series. 

Three other facts stand out. First, the steep inner re- 
gion of the surface brightness profiles is systematically equal 
to or shorter than the isophotally thick part of the galax- 
ies. Second, the isophotally thick part is itself systematically 
contained within the fiat intermediate region of the surface 
brightness profiles. Third, the steep inner region of the sur- 
face brightness profiles is much more prominent along the 
major-axis than in the vertically-summed profiles. 

We are thus led to radically alter the classic 'bulge + 
disc' model, composed of a thick and steep spheroidal bulge 
largely decoupled from a thin (possibly barred) exponential 
disc. Analogously to Athanassoula (20qJ), we propose here 
that galaxies with a B/PS bulge are composed of a thin con- 
centrated disc (a disc- like bulge), formed secularly by the 
bar and responsible for the steep inner region of the surface 
brightness profiles, contained within a (partially) thick bar 
(the B/PS bulge), responsible for the flat intermediate region 
of the surface brightness profiles and the complex morpho- 
logical structures, itself contained within a thin outer ex- 
ponential disc (see Figure]^. Those components are closely 
intertwined dynamically and are largely made of the same 
(disc) material, shaped over long timescales by the bar. 

The challenge to any competing formation scenario for 
galaxies with a B/PS bulge, which represent at least 45 per 
cent of the local galaxy population (iLiitticke et al l2000 a1. IS 
thus to simultaneously and self-consistently explain, equally 
well or better, their numerous morphological, photometric, 
and kinematic properties, as well as the correlations between 
them. 
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